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A scheme for conditional generating a Hermite polynomial excited squeezed vacuum states (HESVS) 
is proposed. Injecting a two-mode squeezed vacuum state (TMSVS) into a beam splitter (BS) and 
counting the photons in one of the output channels, the conditional state in the other output channel 
is just a HESVS. To exhibit a number of nonclassical effects and non-Guassianity, we mainly investigate 
the photon number distribution, sub-Poissonian distribution, quadrature component distribution, and 
quasi-probability distribution of the HPESVS. We find that its nonclassicality closely relates to the con¬ 
trol parameter of the BS, the squeezed parameter of the TMSVS, and the photon number of conditional 
measurement. These further demonstrate that performing the conditional measurement on a BS is an 
effective approach to generate non-Guassian state. 
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1. INTRODUCTION 

Quantum state engineering has been a subject of in¬ 
creasing interest to construct various novel nonclassi¬ 
cal states in quantum optics and quantum information 
processing dll, 0]. From a theoretical point of view, 
the simplest way of generating nonclassical field states 
is to apply the photon creation operation to classi¬ 
cal states such as the thermal and coherent states P3I- 
i]. These nonclassical states, such as the single-photon 
added coherent state ll^ and single-photon-added ther¬ 
mal state 10], have been realized experimentally. Sub¬ 
sequently, it has been demonstrated that subtracting 
photons from traditional quantum states exhibit an 
abundance of nonclassical properties lIsl- flTh . Photon 
subtraction or addition can improve entanglement be¬ 
tween Guassian states ide , loophole-free tests of Bell’s 
inequality idl], and quantum computing ifT^I . 

To meet the requirement of the development of quan¬ 
tum optics and quantum information tasks, some non¬ 
classical states are explored by performing the dif¬ 
ferent combination of photon subtraction and photon 
addition lfl5l - [2Qh . which have different properties. Kim et 
al ll^l discussed single photon adding then subtracting 
(or single photon subtracting then adding) coherent state 
(or thermal state) to probe quantum commutation rules 
[a,a^] = 1. Lee et al flT^ investigated the nonclassical¬ 
ity of field states when photon subtraction-then-addition 
operation or the photon addition-then-subtraction oper¬ 
ation is applied to the coherent state (or thermal state), 
respectively. Yang and Li llT^ analyzed multiphoton ad¬ 
dition followed by multiphoton subtraction (a^a^^) and 
its inverse (a^^a^) on an arbitrary state. Recently, Lee and 


Nha ll2^ proposed a coherent superposition of photon 
addition and subtraction, ta + ra^ (|t|^ + |r|^ = 1) acting 
on a coherent state and a thermal state. More recently, 
we investigated the nonclassical properties of optical 
fields generated by Hermite-excited coherent state ll2^ 
and Hermite-excited squeezed thermal states JH]. 

On the other hand, another promising method for 
generating highly nonclassical states o f op tical fields is 
known to be conditional measurement ll26l43Qh . Namely, 
when a system is prepared in an entangled state of two 
subsystems and a measurement is performed on one sub¬ 
system, then the quantum state of the other subsystem 
can be reduced to a new state. In particular, it turned 
out that conditional measurement on a beam splitter may 
be advantageously used for generating new classes of 
quantum states ll^, 1^1 . Dakna’s group used conditional 
measurement on the BS to generate cat-like state ifsill . 
Podoshvedov et al|[^ proposed optical scheme for gen¬ 
erating both a displaced photon and a displaced qubit 
via conditional measurement. In Ref, ll^ . they proposed 
to create arbitrary Fock states via conditional measure¬ 
ment on the BS. In addition, conditional output measure¬ 
ment on the BS may be used to produce photon-added 
states for a large class of signal-mode quantum states, 
such as thermal state, coherent state, and squeezed 
states (UK]. Similarly, photon-subtracted states can be 
produced by means of conditional measurement on the 
BS lIs^ . Therefore, based on conditional measurement 
on the BS, it is possible to generate and manipulate var¬ 
ious nonclassical optical fields in a real laboratory. 

In this paper, we study the Hermite polynomial excited 
squeezed vacuum state (HESVS), a kind non-Gaussian 
quantum state, generated by conditional output mea- 
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surement on a BS. The calculations show that when a 
two-mode squeezed vacuum state (TMSVS) is injected in 
the input channels and the photon number of the mode 
in one of the output channels is measured, then the mode 
in the other output channel is prepared in a conditional 
state that has the typical features of a Hermite polyno¬ 
mial excited squeezed state. To exhibit the nonclassi- 
cal properties of this conditional state, we mainly ana¬ 
lyze the states in terms of the photon number distribu¬ 
tion, sub-Poissonian distribution, quadrature component 
distribution, and Quasi-probability distribution including 
the Wigner function(WF) and Husimi function(HF). The 
paper is organized as follows. Section 2 presents the ba¬ 
sic scheme for generation of the HESVS and its normal¬ 
ization related to Legendre polynomial. The nonclassical 
properties of the HESVS are analytically and numerically 
studied in Section 3-4. The results indicate that the con¬ 
ditional HESVS is strongly noncassical and non-Gaussian 
due to the presence of the partial negative WF. Finally, a 
summary and concluding remarks are given in Section 5. 



FIG. 1: Preparation scheme of HPESVS. When a TMSVS is 
mixed by a beam splitter and the number of photons |m)is 
measured in one of the output channels, then the conditional 
quantum state in the other output channel is generated. 


II. GENERATION OF HERMITE POLYNOMIAL EXCITED 
SQUEEZED STATE 


operator S 2 (r) on the two-mode vacuum. 


It is well known that the input-output relations at a 
lossless beam splitter can be characterized by the SU(2) 
Lie algebra. In the Schrodinger picture, the role played 
by the beam splitter (BS) upon the input state results 
in the output state 


Pout = ( 1 ) 

where B = exp \0 (a^6 — ab^)] corresponds to the uni¬ 
tary operator in terms of the creation (annihilation) op¬ 
erator (a) and 6^ (b) for mode a and b, whose transfor¬ 
mations satisfying] 


BaW = acosO — 6sin6>, 

BbB^ = a sin 0 b cos 0. (2) 

Moreover, cosO and sinO are the transmittance and re¬ 
flectance of the beam splitter, respectively. Note that the 
globe phase factor of BS may be omitted without loss of 
generality. For the sake of simplicity, we also assume that 
0 is tunable in the range of [0,7r/2]. Under special cir¬ 
cumstances, when 6> = 0or6> = 7r/2, the BS corresponds 
to the cases of total transmission and total reflection, re¬ 
spectively. For 0 = tt/A, the BS is just the symmetrical, 
i.e. 50/50 BS. 


A. Hermite polynomial excited squeezed state 

A two-mode squeezed vacuum state (TMSVS) is the 
correlated state of two field modes a and b (signal and 
idle) that can be generated by a nonlinear medium. The¬ 
oretically, the TMSVS is obtained by applying the unitary 


I^')a6 = ^2 (r) |0,0) = cosh-1 |0,0), (3) 


where £'2 (r) = ex.p [r — ab)] is the two-mode 
squeezed operator and the values of r determines the 
degree of squeezing. The larger r, the more the state is 
squeezed. Especially, when r = 0, reduces to two¬ 
mode vacuum state |0,0). 

The conceptual scheme of the experimental setup is 
depicted in Fig.l. The two input modes prepared in the 
two-mode squeezed state (p^^ = 1^)^^ (^|) is mixed at 
BS, so the output-state density operator can be given by 
Pout = B {4'| St. In fact the output modes in 
are generally highly correlated. When the photon num¬ 
ber of the mode in the second output channel is mea¬ 
sured and m photons are detected, then the mode in 
the first output channel is prepared in a quantum state, 
whose density operator reads as 


Pout = b H B {^1 St \m)^ = , 

( 4 ) 

where is the normalized output conditional state 
(a pure state) and is the normalization factor deter- 
mined by Tr = 1. 

Next, using the integration of the TMSVS lls^i . 


I^)a6 = 


1 f cPa 
sinh r J tt 


|a//tanh r+aa^+a*6^ 


lOajOft) , 


( 5 ) 


where |0a,0b) = |0)^ ® |0)j is two-mode vacuum state, 
and the transformation relation in Eq.©, after some al- 
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gebra we derive that 


b {m\ B 1^') 


ab 


[ 

sinhr J 
1 d™ 


72 

^ Ictp/tanh rg(a cos 0+Ck:* sin 9)a^ |q^ 

TT ^ 

^r{a* cos 9—a sin 0) i 

r=0 


dB 


coshrVml dr^ 


-62 


-\-rua^ 


|0)a k=0, 


( 6 ) 


where fi = sin26>tanhr and v = cos26>tanhr. Hence the 
output conditional state is explicitly expressed as 


5i (A) |0) 


(7) 


with Q.rn = cosh A/(2’^m!A^rn cosh^ r), where we have 
used the generating function of the single-variable m- 
order Hermite pol 3 momial Hm (x) = d!pe^^'^~^^\r=o and 
the expression of single-mode squeezed vacuum 
S'! (A) |0) = cosh“^/^ |0) with the single¬ 

mode squeezed operator Si (A) = exp [A — a^) /2]. 
Eq.([7]) indicates that the conditional state j^^n) is actu¬ 
ally a single-mode m-order Hermite polynomial excited 
squeezed vacuum state. It is worth noticing that the de¬ 
gree of squeezing A of the conditional state is not the 
same squeezing parameter r of the TMSVS but related 
to the parameters of the TMSVS and the BS satisfying 
tanhA = sin26>tanhr. For the symmetrical case, i.e., 
e = B = Sia (r) |0) ® Sib (-r) |0) is just 

the product state of two separate single-mode squeezed 
vacuum state. In this case, when we detect m photons 
in the second output channel, the conditional state is al¬ 
ways Si (r) |0) with the same squeezed parameter r. This 
is that the effect of the symmetrical BS splits the entan¬ 
gled TMSVS into two independent single-mode squeezed 
vacuum state. Note that when no photons are detected, 
m = 0, then also reduces to Si (A) |0). 


B. Normalization via probability of such event 


In addition, the normalization factor Nm is determined 
by the probability p (m) of such an event given by 


Nm = p{m) 

= Tr(^b{m\B\^)^,{^\BUm),) 

Bis^l2-BiT^I2+B2ST 


d2r, 


mlv^cosh^ r ds^^dr^ 
cosh^ r^/A 


{•/Bl) 


s=r=0 


( 8 ) 


where we have set ^4 = I — jj?, Bi = 11 / (Acosher), 
B 2 = v'^IA, B 3 = (tanhV - Ai") jA, B^ = j (A^bJ), 



FIG. 2: The probability of producing |^rn) is shown as a func¬ 
tion of the parameter r of the TMSVS for two parameter values 
of the BS [(a) 0 = tt/T; (b) 0 = 27r/7] and various values of 
m, where m = 1,2,3,4correspond to the solid, dashed, dotted 
and dotdashed lines, respectively. 


and in the last step we have used the formula of m-order 
Legendre polynomial (x), i.e., 

2™to! „ , , 

dt^dT^^ " ' 

(9) 

Especially when r = 0, \^) |0,0), there is no photons 

in the output channels. In this case, there is no neces¬ 
sary to making conditional measurement. So the event 
is happen only for r / 0. When = 0 or 6> = 7r/2, 
leading to p = 0 and u = tanh r or u = — tanh r then 
^4 = 1, = 0, ^2 = tanh^ r, Bs = tanh^ r, and = 1, 

Nm\e=Oor 0 = 7 v /2 = tauh^™ r/cosh^ r and |4'm) is just the 
Fock |m), which is rational because of the inherent prop¬ 
erties of the TMSVS. If 6> = 7r/4, the BS is just the sym¬ 
metrical, i.e., 50/50 BS, leading to p = tanhr, u = 0 then 
A = 1 - tanh^ r, Bi = ^ output 

states is the the product of two independent single-mode 
SVS. 

According to Eq.®, we discuss the probability of ob¬ 
serving such a conditional m-order HPESVS. In Fig.2 the 
probability p (m) is plotted for two parameter values of 
the BS. For a given transmittance of BS, p (m) as a func¬ 
tion of the input squeezing parameter r can attain a max¬ 
imum and the maximum is shifted towards larger values 
of r when m is increased (see Fig.2a). In this figure, 
for each r, we should use the transmittance of the BS 
in a way to optimize the success probability. The ideal 
procedure would use the appropriate transmittance for 
each value of r and m. By tuning the parameters of 
the interaction, namely, the control parameter of the BS, 
the squeezed parameter r of the TMSVS, and the pho¬ 
ton number of conditional measurement m, the HPESVS 
may be modulated, generating a wide range of nonclas- 
sical phenomena, as described below. 


III. OBSERVABLE NONCLASSICAL EFFECTS OF THE 
CONDITIONAL HPESVS 

To study the nonclassical properties of the conditional 
states in more detail, we shall calculate the photon 
number distribution, sub-Poissonian distribution, and 
quadrature component distribution. 
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FIG. 3: Photon-number distribution of the conditional HPESVS 
for (a) m = 1,6> = 27r/7, r = 0.5; (b) m = 1, 6> = 27r/7, r = 1.0; 
(c) m = 4, ^ = 27Tl7,r = 1.0; (d) m — jl — 1.0, 

respectively. 


A. Photon number distribution 

The photon number distribution (PND), the probabil¬ 
ity of finding n photons, is a key characteristic of every 
quantum state. Recalling Eq.([7|), the PND of the condi¬ 
tional HPESVS reads as 


P(n|m) = |(n|4'„)|\ (10) 

Using the unnormalized coherent \z) = exp[ 2 ;a^] |0), 
leading to \n) = k) \z=o, and combining with 

Eq.O, we finally obtain 


P {n\m) 


1 


mlnlNm cosh r 
mini 


Q2n 


dr'^ds'^ 




Nm cosh^ r 


mm m,n 


E 

^=0 


„ I m — g m-\-n — 2g 

’ ' (-1)~ (f) " 




, ( 11 ) 


where we have used = m\5mn (•^mn is 

Krocher function) and is given in Eq.®. In the sum¬ 
mation of Eq. ifTTl) . the value of g must make and 
be integer. To see clearly the variation of the PND, 
in Fig. 3 we plot the bar graph of the PND for the condi¬ 
tional HPESVS with different values of parameters m, 0, 
and r. From Fig.3 we easily see that when the number 
m of the conditional measurement is odd (even), then 
the photon-number distribution is nonzero only for odd 
(even) photon numbers. The probability P {n\m) with 
different parity between m and n is zero. For given m 
and 0, the bigger the squeezing parameter r, the wider 
the distribution (see Figs.3(a) and 3(b)). 


B. Sub-Poissonian distribution 


In order to study the photon-number statistics of 
this conditional state, we first calculate = 

(T^rnl Using the completeness of coherent 

state / ^ \a) {a\ = 1 as well as = ^e^^|t=o yields 


1 f)2m 

^ _ t _ ^ ^-Bis^/2-Bit^/2+B2TS 

Nmrnl cosh^ r^/A ds^dr^ 

^ ^ ^ ^[tiiy(xs-\-yr)-\-iy(ys-\-xr)-\-ti(x‘^-\-y‘^)/2-\-xy]/A\_ 

dx^dy^ |s-r-ai-i/- 

( 12 ) 


Thus, the mean photon number 


(n) = (a^a) = (aa^) — 1, (13) 

can be determined by Eq. ([T2l) with k = I = 1. Exam¬ 
ples are shown in Figs.4(a) and 4(c). We see that when 
0 = 7TIb, the number of photons that can be found in 
increases with m for the given larger r, and (n) 
as a function of 0 is symmetric distribution for 0 = tt/A. 
This is simply a consequence of the BS transformation. 
In particular, when no photons are detected, m = 0, then 
(n) = sinh^ r reduces to the mean photon number of 
single-mode squeezed vacuum state. 

A measure of the deviation of the photon number dis¬ 
tribution from a Poissonian is the Mandel Q factor de¬ 
fined bv ll^ 


Q = 


(n2> - {nf 
in) 

— 2 (aa'*') + 1 
(aat) — 1 


(14) 


It holds that Q ^ 0 and the equality is achieved for the 
Fock state. The light is sub-Poissonian when the photon- 
number variance (n^) — (n)^ is less than (n). This is 
indicated by a negative value of Q. The statistics are 
Poissonian when Q = 0, and super- (sub-) Poissonian if 
Q > 0 (Q < 0). 

According to Eqs. ([T^ and ([T4l) , we plot the variation 
of Q for HPESVS versus r or 0 for different m = 1,2,3,4 
in Fig.4(b) and 4(d). It is clearly seen that Q as a func¬ 
tion of 0 is also symmetric distribution for 0 = tt/A and 
the HPESVS has sub-Poissonian statistics behavior due to 
the emergence of the negativity of Q. With the increasing 
values of m, the increasing intrend of Q is accelerated for 
larger r. To further exhibit the high nonclassicality. Fig. 5 
shows the dependence of this conditional state on 0 and 
r for four different Q factors. Especially, we consider first 
the boundary case of the Poissonian distribution, Q = 0 
(see the solid line in Fig.5). 
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FIG. 4: (a) Mean photon number (n) and (b) Mandel Q factor 
versus r for different m = 1 (solid line), m = 2(dashed line), 
m = 3(dotted line), and m = 4(dotdashed line) with the same 
^ = 7r/5. (c) Mean photon number (n) and (d) Mandel Q factor 
versus 0 for different m = 1 (solid line), m = 2(dashed line), 
m = 3(dotted line), and m = 4(dotdashed line) with the same 
r = 0.5. 



C. Quadrature component distribution 

Next we pay attention to the conditional quadrature 
component distribution (QCD) ll^ 

P = \{x,(f \^m)f , (15) 

which can be measured in balanced homodyne detection. 
Here \x, ip) is the eigenstate of the quadrature component 
X {ip) = expressed as in the Fock basis 

|a;, + |g^ _ Qg) 

Using Eqs.O and ([T6l) and inserting the completeness 
of coherent state / ^ |<^) (<^| = after integration, the 
wave function (x, ip reads 


FIG. 5: The sub-Poissonian properties of the conditional state 
with Q = 0(solid line), Q = —0.2(dashed line), Q = 
—0.5(dotted line), and Q = — 0.8(dotdashed line) in the plane 
space of two parameters (^and r) with different number of con¬ 
ditional measurement: (a) m = 1; (b) m = 2; (c) m = 3; (d) 
m = 4, respectively. 


IV. QUASI-PROBABILITY DISTRIBUTION OF THE 
CONDITIONAL HPESVS 

Quasi-probability distribution function in the phase 
space is a very useful tool for a comprehensive descrip¬ 
tion of the nonclassical state. Thus, in this section, we 
shall analytically discuss several quasi-probability dis¬ 
tributions, including Wigner function and Husimi func¬ 
tion to characterize the nonclassicality of the conditional 
HPESVS. 


A. Wigner function 


,-'/*(vr72)"e-»-’ /A \ 

(1 + coshr-^” I XtV ’ 

(17) 

where we have set 0 = l + Tt = (l - /0, 

r = (/i + tanh^ r) /0, and A = e-^^z//0. 

As a result of Eq.([T7l), we easily obtain the condi¬ 
tional QCD defined by Eq. ([T5l) and plot the variation of 
P {x, ip\m) for the HPESVS as a function oix or ip for dif¬ 
ferent m = 1,2,3,4 in Fig.6. Ones see that for ip near 7r/2 
the QCD P (x, ip\m) with m > 0 exhibits two separated 
peaks, wheresas for ip close to 0 and tt an interference 
pattern is observed. 


The WE was first introduced by Wigner in 1932 to 
calculate quantum correction to a classical distribu¬ 
tion function of a quantum-mechanical system. The 
presence of negativity of the WE is a signature of its 
nonclassicality ijHI, Is^ . For a single-mode density op¬ 
erator p, the WE in the coherent state representation \z) 
can be expressed as 



where a = {x ip) / a/2. The Wigner function of the 
conditional state be calculated in 
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FIG. 6: Quadrature-component distribution P(x,(f\m)oi the 
conditional state |^m)for 0 = 7r/7,r = 0.5and various num¬ 
bers mof measured photons with (a) m = 1, (b) m = 2, (c) 
m = 3, and (d) m = 4, respectively. 





FIG. 7: Wigner functions W(x,p\m)of the conditional state 
I^rn) for 0 = 7r/7,r = 0.5and various numbers m of measured 
photons with (a) m = 1, (b) m = 2, (c) m = 3, (d) m = 4, 
respectively. 


a straightforward way. 


B. Husimi function 


W{x,p\m) 




TrNmTnl cosh^ r\/]4 

,, ^Rs+R*r-B2rs-^s^-^P\ 

X ^ _ . e 2 2 |s=T- = 0 


dr^ds^ 
2m! 


o1*^1^ I P ^ I 


ttA^^ cosh^ r^/A 


E 

z=o 


(-^2)^ (-El/2) 
/! [(m — /)!]^ 


m—l 


Hm-l I - 




(19) 


where we have set E = (l + /i^) / (l — /i^) and 
R = 2v (oL - [loP) /A. Especially, when no pho¬ 
tons are detected, m = 0, then W{x,p\0) 

exp (—/tt is a Gaussian form in phase 
space, which is just WF of single-mode SVS, as expected. 


The WFs of the conditional HPESVS in Fig. 7 are 
plotted for the same parameters as in Fig.6. The figures 
indicate that the conditional HPESVS is a noncassical 
non-Gaussian state, since the partial negative regions in 
phase space are observed in Fig. 7. This further demon¬ 
strates that performing the conditional output measure¬ 
ment on a BS is an effective approach to generate non- 
Guassian state. In addition, it is seen from Fig.7 that for 
odd m there exists a negative valley in the center region, 
whereas for even m there exists a main peak. In fact, for 
the center region 147(0,0|m) = | (-1)’^, as expected. 


The Husimi function Q{x,p\m) of the state is de¬ 
fined bv ll^ 


Q{x,p\m) =-\{/3 \'^rn)\^ , (20) 

TT 

where |/3) is a coherent state and (3 = {x ip) /\f2. Us¬ 
ing Eq.©, the scalar product (/3 can be easily cal¬ 
culated as follow 


(/^ \^m) 




on 




yJmlNjn coshr 


r=0 1 


( 21 ) 


So we find that Q {x,p\m) can be written as 


Q{x,p\m) 


TTNrnm}. cosh^ r 


Hr, 



( 22 ) 


As expected, for m = 0 the Husimi function is Gaus¬ 
sian, while for odd m a two-peak structure and for even 
m a single peak are observed in Fig.8, whose parameters 
are the same as WFs in Fig.7. Note that the Husimi func¬ 
tion is a phase-space function that can be measured in 
multiport balanced homodyning. Since the Husimi func¬ 
tion can be regarded as a smoothed Wigner function, it is 
always non-negative and the oscillating behavior, typical 
of WF (see Fig.7), cannot be observed. 
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ditional measurements using a simple beam splitter 
scheme. When a two-mode squeezed vacuum state is 
mixed by a beam splitter and the number of photons 
is measured in one of the output channels, then the 
conditional quantum state in the other output channel 
reveals all properties of a Hermite polynomial excited 
squeezed state. Then, we also have numerically analyzed 
the conditional HPESVS in terms of the photon-number 
statistics, quadrature-component distribution and quasi¬ 
probability distribution such as the Wigner and Husimi 
functions. The results show that by tuning the param¬ 
eters of the interaction, namely, the control parameter 
of the BS, the squeezed parameter r of the TMSVS, and 
the photon number of conditional measurement m, the 
HPESVS may be modulated, generating a wide range of 
nonclassical phenomena. This further demonstrates that 
performing the conditional measurement on a BS is an 
effective approach to generate non-Guassian state. 


FIG. 8: Husimi functions Q(x,p\m) of the conditional state 
l^m) for 0=7rl7, r = 0.5and various numbers m of measured 
photons with (a)m = 1, (b) m = 2, (c) m = 3, (d) m = 4, re¬ 
spectively. 
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